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ABSTRACT 

We present direct upper limits on gravitational wave emission from the Crab pulsar using data from 
the first nine months of the fifth science run of the Laser Interferometer Gravitational- wave Observa- 
tory (LIGO). These limits are based on two searches. In the first we assume that the gravitational 
wave emission follows the observed radio timing, giving an upper limit on gravitational wave emission 
that beats indirect limits inferred from the spin-down and braking index of the pulsar and the ener- 
getics of the nebula. In the second we allow for a small mismatch between the gravitational and radio 
signal frequencies and interpret our results in the context of two possible gravitational wave emission 
mechanisms. 

Subject headings: gravitational waves - pulsars: individual (Crab pulsar) 
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1. INTRODUCTION 

The Crab pulsar (PSRB0531+21, PSR J0534+2200) 
has long been regarded as one of the most promis- 
ing known local sources of gravitational wave emission 
and is an ic onic target for gravitatio nal wa ve searches 
(Press fc Thorne 1972; Zimmermani] I1978D . Its high 
spin-down rate, v ~ —3.7x10^^° Hz s~^, corresponds to a 
kinetic energy loss rate oi E = 4:7r^Izzi^\i^\ ~ 4.4x10^^ W 
(using a spin frequency of = 29.78 Hz and the canonical 
value of lO^^kgm^ for the principal moment of inertia 
Izz-) This loss is due to a variety of mechanisms, in- 
cluding magnetic dipole radiation, particle acceleration 
in the magnetosphere, and gravitational radiation. If 
one assumes that all the energy is being radiated grav- 
itationally, the gravitational wave tensor amplitude at 
Earth is hf = SMxlQ-^^ l3sr^^^{\i>\/iy)^/'^, where rkpc 
is the distance to the pulsar in kpc and I^s is the moment 
of ine rta in units of the canonical value (jAbbott et al.l 
I2007c[ ). For the Crab pulsar this "spin-down upper limit" 
is hl*^ — 1.4 X 10~^^, using the canonical moment of in- 
ertia and a distance r = 2 kpc. It has long been known 
that the Laser Interferometer Gravitational-wave Obser- 
vatory (LIGO) can achieve this sensitivity by integrating 
several months of data with the initial design noise spec- 
trum. 

The electromagnetic emission and accelerating expan- 
sion of the Crab Nebula are powered almost entirely by 
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the rotation of the pulsar. The question now is whether 
these two loss mechanisms can account for the vast ma- 
jority of the observed rotational energy loss, or whether 
gravitational wave emission has a significant part to play. 

The bolometric luminosity of the nebula is (1- 
2)xl0'^^ W, which accounts for less than half the spin- 
down power fe.g.. [Davidson fc FesenI ITOSSl ) . There have 
been many attempts to estimate the power involved in 
the observed a cceleration of optical filaments , for exam- 
ple recently bv lBeiger fc Haensell (|200H 2003'). However 
these depend on poorly known factors such as the mass 
and expansion history of the nebula, and the uncertain- 
ties in the estimated power are comparable to the spin- 
down power itself. Thus electromagnetic observations of 
the nebula, within their uncertainties, still allow for a 
substantial fraction of the spin-down power to be emit- 
ted in gravitational waves. 

The braking index n = vi) /v^ of the pulsar further con- 
strains the gravitational wave emission. The observed 
value n — 2.5 still is not well understood on theoretical 
grounds, but since quadrupolar radiation has n = 5 it im- 
plies that only a small fraction of the spin-down power 
is emitted in gravitationa l waves. The best estimate in 
print is bv lPalombal (|2000D . who uses a phenomenological 
model of the spin-down (present and historical) together 
with the present braking index and known age of the 
pulsar to estimate that the highest possible ho today is 
about 40% of the spin-down limit. This value is con- 
sistent with the observations of the nebula, and is also 
observable with several months of data from LIGO's fifth 
science run (S5). 

Early directed searches for gravit ational waves from 
the C rab pulsar were performed by iLevine fc Stebbini 
(|1972t ). using a 30 m i nterferometer to give a strain upper 
limit of 3xl0~^^, and lHirakawa et all ('1978'), using a bar 
detector. The most recent bar result (Suzuki 1995) gave 
an upper limit that was still over an order of magnitude 
above the spin-down limit. The LIGO detectors have 
improved on these results, with LIGO's second science 
run (S 2) producing a 95% upper limit of h^^^" = 4.1 x 
10"^^ () Abbott et al.ll2005f) . and the combined data from 
the S3 and S4 runs produced an upper hmit of h^^'^" — 
3.1 X 10-24 ([Abbott et al.ll2007cD only 2.2 times greater 
than the spin-down limit. 

In this Letter, we describe searches of data from the 
fifth LIGO science run, which star ted on 2005 Novembe r 
4 and ended on 2007 October 1 (| Abbott et al.ll2007bD . 
During this period the detectors (the 4 km and 2 km de- 
tectors at LIGO Hanford Observatory, HI and H2, and 
the 4 km detector at the LIGO Livingston Observatory, 
LI) were at their design sensitivities and had duty fac- 
tors of 78% for HI 79% for H2, and ~ 66% for LI. The 
GEO600 detector (jLiick et al.ll2006[ ) also participated in 
the S5 run but was much less sensitive at the frequency 
of the expected signal. 

The Crab pulsar was observed t o glitch on 2006 August 
23 at approximately 04:00 UTC (jLvne et al.ll2007l: iLvnd 
[2006h . Since the glitch mechanism is not certain and 
may involve unpredictable changes in the gravitational 
wave timing and amplitude, we use this glitch as natural 
point at which to pause this coherent search for the Crab 
pulsar. Our data set consists of HI and H2 data from 
2005 November 4 and LI data from 2005 November 14 up 



to 2006 August 23. For the two different searches carried 
out in this analysis, described below, this gives 201, 222 
and 158 days of data for HI, H2 and LI respectively 
for the single-template search, and 182, 206, and 141 
days of data for the multi-template frequency-frequency 
first derivative search, which required larger contiguous 
segments than the single-template search. 

2. METHODS 

We use two different methods (see I Abbott et al1l2004j ) 
to search for gravitational waves from the Crab pulsar 
to account for different emission scenarios. One method 
uses a single time domain template for the gravitational 
wave signal assuming that the gravitational wave pe- 
riod evolves precisely as the electromagnetic pulse pe- 
riod. The other method works in the frequency domain 
to cover a relatively small, physically motivated range 
of frequency and spin-down values. The searches use 
the known frequency and position of the Crab pulsar, 
as derived from the Jo drcU Bank Crab Pulsar Monthly 
Ephemeris (|Lvne et al.ir20 07). Using this ephemeris and 
the assumption that the gravitational wave and electro- 
magnetic phase track each other precisely, we can predict 
the signal phase evolution with negligible uncertainty. 
Both searches assume that emission will be at or near 
twice the pulsar's spin frequency, 2v — i^gw ~ 59.56 Hz, 
which is the frequency of emission by a steadily rotating 
quadrupolar deformation, i.e. a triaxial star. The Crab 
pulsar might be emittin g at z/qw ~ 4z//3 through an r- 
modc (jOwen et al.lll998[ ) if the mode saturates at a small 
amplitude and thus is long-lived (e.g.. Brink e t al« 2005) . 
However, the uncertainty of this frequency is relativel y 
large, of order one part in 10^ (jLindblom et al.iri999f) . 
Due to this, and the greater instrument noise at this fre- 
quency, we did not search for r-modes. Although 2iy is 
close to the 60 Hz power line frequency, it is sufficiently 
far away that the searches are relatively unaffected by 
non-stationary components of the power line noise. The 
absolute timing accuracy of the LIGO data is sufficiently 
good that the likelihoods produced for each detector can 
be combined to give a joint likelihood. 

For a given search frequency and spin-down, the 
four unknown signal parameters are the gravitational 
wave amplitude ho, the initial phase (j)o, the spin-axis 
inclination angle l, and the polarization angle ip. X-ray 
observations of the Crab Pulsar Wind Nebula provide 
values of the orientat ion angle l and polarisation angle ^p 
of the pulsar. From iNg fc Romaiiil (|2004l . [200l) we use 
t = 62.17 ± 2.195° and V = 125.155 ± 1.355°, where we 
have taken the mean of the best fit values for the outer 
and inner tori of the nebula. We use these ranges to 
put Gaussian priors on these two parameters for both 
the search techniques. On the chance that the star is 
misaligned from these structures, we also present results 
using uniform priors over the allowed ranges of the pa- 
rameters. 

The single-template search (jPupuis fc Woan|[2005l) as- 
sumes a triaxial star emitting gravitational waves at 
precisely twice the spin frequency, following the elec- 
tromagnetic pulse phase evolution and taking into ac- 
count the small variations in phase caused by timing 
noise (|Pitkin fc Woan 20071. It uses a standard Bayesian 
methodology to produce a joint posterior probability vol- 
ume over the four unknown parameters using data from 
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all three detectors. We use both uniform priors and re- 
stricted priors on and t when calculating the posterior. 
We marginalize the angle parameters to produce a pos- 
terior probability for ho and from this calculate a 95% 
degree-of-belief upper limit on the gravitational wave am- 
plitude. 

A search was also performed at gravitational wave fre- 
quencies z^GW in a narrow band about 2v, based on 
simple astrophysical arguments. We begin by writing 
'^GW = 2i^(l + S), where 5 is a small number. A rela- 
tion of this form holds if the gravitational waves are pro- 
duced by a component spinning separately from the elec- 
tromagnetically emitting one, with the two components 
linked by some torque which acts to enforce co-rotation 
between them on a timescale TcoupUng- In such a case 

^ ^ ^coupling/ '^spin— down: 

where Tspin-down ~ I'/i' - 2500 
years. A relation of the form given for i^gw above also 
holds if the gravitational waves ar e produced by free 
prece ssion of a nearly biaxial star (jJones fc AnderssonI 
I2OO2I ). In such a case 5 ^ a{Izz — Ixx)llxx where a is a 
factor of order unity dependent on the geometry of the 
free precession, e.g. the angle between the symmetry axis 
and angular momentum axis. No clear signature of free 
precession has been seen in the radio pulsations of the 
Crab pulsar, although precession would have little effect 
on the radio signal if the amplitude of the precession were 
small. 

Together, these scenarios suggest searching over a fre- 
quency interval iAi^GW centred on where Ai/GW 
\5\ 2v. We have followed such a strategy, using a 
maximum value of |^| = 10"'*. In terms of the 
two-component model, such a |^| value corresponds to 
7"coupiing ~ 10""^ Tspin-down ~ scvcral months, compara- 
ble to the longest timescales seen in glitch recovery where 
re-coupling between the two components might be ex- 
pected to occur. In terms of free precession, \5\ = 1G~^ 
is on the high end of the range of deformations that 
com pact objects are thought to be capable o f sustain- 
ing JOwed[200l[Li^[2007l : lHaskell et al.l[200l . 

Using the above estimates as a guide, a band of fre- 
quencies ±6 X 10^^ Hz centred on twice the Crab pulsar's 
observed frequency was searched over. Corresponding 
bands in frequency derivatives were motivated via differ- 
entiation of the equation for i^gw, which together with 
the assumption that S itself evolves no more rapidly than 
on the spin-down timescale, leads to a band in frequency 
first derivative of ±1.5 x 10"^'^ Hz/s, with searches over 
higher derivatives being unnecessary. 

The multi-template search method is a maximum like- 
lihood tec hnique, the coherent mul ti-detector F-statistic 
derived in ICutler fc Schut j (|2005( ) . An explicit search 
is required over a single sky position and second deriva- 
tive of the frequency, and over the selected ranges of the 
frequency and of the first frequency derivative. The spac- 
ing of the templates is chosen in such a way as to ensure 
at most a 5% loss in the detection statistic, resulting 
in a total of 3 x 10^ templates. The detection statis- 
tic 2F is computed for each template. The expected 3a 
range of the largest 2F value for Gaussian noise (no sig- 
nal present) and 3 x 10^ templates is 35-49. The largest 
2F value found in the actual search is 37, well within the 
expected range for noise. 

Based on the largest 2F value, 95% confidence upper 
limits are produced using a frequcntist Monte Carlo in- 



jection method, as described in lAbbott et al.l ()2007af) . 
For the unknown parameters uniform distributions and 
physically informed distributions were used for the in- 
jected population of signals, consistent with the choices 
made for the single-template time domain search. 

3. RESULTS 

In the single-template search the joint (i.e. multi- 
detector) posterior probability distribution for the grav- 
itational wave amplitude peaks at zero, indicating that 
no signal is visible at our current sensitivity. The joint 
95% upper limit on the gravitational wave amplitude, 
using uniform priors on all the parameters, is /iq^'^° = 
3.4 X 10~^^. In terms of th e pulsar's ellipticity, g iven by 
e = 0.237 /i_24r-kpcJ^"^/38 (|Abbott et al.l l2007c( l. where 
ft,_24 is Hq in units of IxlO"^**, this gives e = 1.8xl0~* us- 
ing the canonical moment of inertia and r = 2 kpc. This 
is 4.1 times lower than the spin-down upper l imit and 
also 1 .6 times lower than the limit estimated by Palombal 
(|2000D (see gH) Squaring the ratio of the spin-down and 
direct upper limit shows that less than w 6% of the to- 
tal power available from spin-down is being emitted as 
gravitational waves, assuming the canonical moment of 
inertia. Using the restricted priors on ip and t wc get 
an upper limit on ho of 2.7 x 10^^''', which is 1.3 times 
smaller than that with uniform priors, and corresponds 
to less than 4% of the spin-down energy available. 

With the coherent multi-template frequency-frequency 
first derivative search we set 95% confidence upper limits 
on ho and ellipticity of 1.7x10"^** and 9.0x10"'* respec- 
tively, over the entire parameter space searched. These 
upper limits are larger than the single-template search 
limits by roughly a factor of five. This is to be expected 
because the larger number of templates raises the num- 
ber of trials and thus the statistical confidence threshold. 
Assuming restricted priors on -0 and t yields an improved 
upper limit of 1.2x10"^'', a factor of 1.2 below the spin- 
down limit, across the entire parameter space searched. 
This limits the energy budget of gravitational waves to 
be less than 73% of the available energy. These quoted 
upper limits are subject to uncertainty in the calibra- 
tion of the detectors. Amplitude calibration uncertain- 
ties for HI, H2 and LI, respectively, are: 8.1%, 7.2% 
and 6.0% (single-template analysis), and 9.5%, 7.8% and 
8.7% (multi-template analysis). 

4. DISCUSSION 

Under the assumption that the gravitational wave and 
the electromagnetic signals are phase-locked, our single- 
template search results constrain the gravitational wave 
luminosity to be less than 6% of the observed spin-down 
luminosity. This beats the indirect limits inferred from 
all electromagnetic observations of the Crab pulsar and 
nebula. 

Our upper limits are interesting because they have en- 
tered the outskirts of the range of theoretical predictions. 
Normal neutron stars are believed to be mostly fluid 
with maximum elastic deformations orders of magnitude 
smaller than the few x 10^* of our upper limits, but some 
theories of quark matter predict solid or mostly solid - 
stars which could sustain su ch eUipticities (lOwenI [20051 : 
lLh][2007l : lHaskell et al.l l2007). However, our upper limits 
do not constrain the composition of the star and cannot 
constrain any fundamental properties of quark matter. 



Beating the Crab pulsar spin-down Umit 




ellipticity 



Fig. 1. — The single template search upper limits from S5, for 
the uniform and restricted prior ranges, and spin-down upper limit 
plotted as exclusion regions in a moment of inertia-ellipticity plane. 
Areas to the right of the diagonal lines are excluded. The dashed 
horizontal lines represent estimates of the theoretical lower and up- 
per bounds of acceptable moments of inertia at (l-3)xl0''^ kgm^. 
The shaded area represents the region that is newly excluded with 
these results. 

The ehipticity is proportional to the quadrupolar strain, 
which may simply be very low for a given star no matter 
its composition. The Crab is likely to have an ellipticity 
at least about 10~^^ due to the stresses of its internal 
magnetic field (|Cutlei| [2002[) if the internal field is com- 
parable to the external dipole of 4 x 10^^ G. Our upper 
limits can be interpreted as direct upper limits of about 
10^^ G on the internal magnetic field, de pending on the 
ratio of toroidal to poloidal components (jColaiuda et al.l 
11)08). 

As discussed in lAbbott et al.l ()2007d l there is consid- 
erable uncertainty in the true value of the Crab pulsar's 
moment of inertia. The best guesses at its value come 
from neutron star equation of state models rather than 
direct measurements. Previous pulsar ellipticity upper 
limits and spin-down limits have made use of the canoni- 
cal value of Izz ■ We can however cast our upper limit in a 
way that makes no assumptions about the moment of in- 
ertia, by placing the limit on the neutron star quadrupole 
moment « Izz^- This then allows us to plot the single- 
template search results as exclusion regions in the I-e 
plane. The results, with uniform and restricted prior 
ranges, are plotted in this way in Figure [T] Our up- 
per limits are smaller than the spin-down limit by a fac- 

1 /2 

tor that varies as Izz ■ If we take the theoretical upper 
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bound on the momen t of inertia to be 3x10'^^ kgm^ as in 
(jAbbott et al.ll2007ch then the result with uniform priors 
beats the spin-down limit by a factor of 7.2. 

Finally, the physical interpretation of our multi- 
template search depends upon the assumed cause of 
the splitting i^gw = 2^(1 -I- 5) between gravitational 
and electromagnetic signals. In the context of the two- 
component spin-down model, our results show that a 
gravitational wave emitting component of the star cou- 
pled to the electromagnetic (radio) emitting component 
on a timescale of a few months or less has a quadrupole 
asymmetry lyy — I^x of no more than 9.0 x 10^'*kgm^. 
This is about five times larger than the bound on lyy—Ixx 
obtained in the single-template search. If free precession 
is responsible for the frequency splitting our results in- 
stead give an upper limit on the product AI sin" 6, where 
A J is the Izz — I XX part of the quadrupole moment ten- 
sor th at participates in the prec ession and 6 the wobble 
ang le jj ones &: Andersso^l2002l ). 
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